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Correlation measurements on the states of two-level atoms having passed through a micromaser at different times can be 
used to infer properties of the quantum state of the radiation field in the cavity. Long(short) correlation length in time is to some 
extent associated with super(sub)-Poissonian photon statistics. The correlation length is also an indicator of a phase structure 
much richer than what is revealed by the usual single-time observables, like the atomic inversion or the Mandel quality factor. 
In realistic experimental situations the correlations may extend over many times the decay time of the cavity. Our assertions are 
verified by comparing theoretical calculations with a high-precision Monte-Carlo simulation of the micromaser system. 
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The one-photon atomic transition micromaser pro- 
vides an impressive experimental realization of the interac- 
tion between a two-level atom and a second quantized single- 
mode electromagnetic field. The microlaser [|J is the counter- 
part in the optical regime. Here we focus our interest on the 
micromaser system for which quantum collapse and signs of 
quantum revival have been observed Q. A superconducting 
niobium cavity, cooled down to a temperature of T = 0.5 K 
(corresponding to a thermal photon occupation number of 
rib = 0.15 at the maser frequency of 21.5 GHz), has been 
used to study the quantum state of the radiation field [|7]j8]]. 
The high quality factor of the cavity corresponds to a photon 
lifetime of T cav = 0.2 s. 

In this context, a basic physical feature is the close con- 
nection between the cavity's steady-state photon statistics and 
the fluctuations in the number of atoms in the lower maser 
level, for a fixed cavity transit time r of the atoms [p 10 1 . 
The experimental results of [^] are clearly consistent with the 
appearance of non-classical, sub-Poissonian statistics of the 
radiation field, and exhibit the intricate relation between the 
atomic beam and the quantum state of the cavity. 

In this letter we study in more detail this relation, and dis- 
cuss the role of correlations for revealing the quantum state 
of the micromaser system. The sequence of outgoing atoms 
is viewed as a one-dimensional binary spin chain, each spin 
representing the state of an atom after the interaction with the 
cavity. The spin average is closely related to what is usually 
called the atomic inversion. From the second-order correla- 
tion functions of the spin chain we are immediately led to the 
physical concept of an atomic correlation length £4, charac- 
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terizing the long-range correlations in time in the outgoing 
atomic beam. 

By means of Monte Carlo simulations we study the dynam- 
ical generation of such correlations in the atomic beam and 
how these would show up in an actual experimental situation. 
We discuss the set of characteristic dynamic phases revealed 
by £4, and the implications for the underlying photon distri- 
bution in the cavity. 

The theory of the micromaser has been developed in [S], 
and we follow the notation of that paper. In the experi- 
ments, excited atoms are injected into the cavity at a rate 
R high enough to pump up the cavity from vacuum, i.e. 
R > 7 = 1 /Tcav, or N > 1 in terms of the dimensionless 
flux variable N = R/j. The atom-field interaction time r 
is so much shorter than the average time between two atoms, 
T = 1/R, that with a very high probability only one atom at 
a time is found inside the cavity fll~i|]. Since r is also much 
shorter than T cav , damping effects may be ignored while the 
atom passes through the cavity. If an atom enters at time t 
and interacts for a time r, then the statistical operator of the 
whole system changes to p{t + r) = e~ T p(t)e %Hr , where 
H is the total Hamiltonian of the atom-field interaction. It 
can be approximated by the Jaynes-Cummings (JC) Hamil- 
tonian [|l2p. This Hamiltonian has the property that it only 
affects the reduced density operator pFif] = Tta (/»(*)) of 
the radiation field along its diagonals, so that if pp (t) is diag- 
onal — which we assume — then so is pp(t + t). Writing 
Pf(£) = J2^LoP n (t)\ n )( n \ we may express the above inter- 
action as a linear operator acting on the infinite-dimensional 
vectorp(f) = {po(t),pi(t), . . .}, i.e. asp(t + r) = M(r)p(t), 
where M — M + + M~ is composed of two parts, represent- 
ing either that the photon occupation number is unchanged 
(the atom is in state +), or that it is increased by 1 due to the 
decay of the atom (the atom is in state — ). From the JC-model 
we have 
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M+ m {r) = {I - q n+1 {T))5 n 
M nm( T ) = Qn(r)S n -l, m , 
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(1) 



where at resonance between the cavity mode and the atomic 
transition q n (r) = sin (fir\/n). The quantity f2 is the single- 
photon Rabi frequency [ p"2| ] . 

Let the next atom arrive in the cavity after a time T 3> r. 
During this interval the cavity damping is described by a con- 
ventional master equation, which also preserves the diagonal 
form of the cavity density matrix. It may be brought to the 
form P (t) — —~fLp(t), where 7 is the characteristic damping 
constant of the cavity and L is the matrix: 

Lnm = {n b + l)(nS nim - (n + l)5 n+ x, m ) 

+ n b ((n + l)6 ntm - nS n -i >m ) . (2) 

The statistical state of the cavity when the next atom arrives is 
thus given by 



p(t + T)= e-~< LT M(T)p(t) 



(3) 



The time intervals between atoms are assumed to be Poisson- 
distributed dV(T) = cxp (-RT)RdT with an average inter- 
val T = 1/R (for a discussion of non-Poissonian beam statis- 



tics see [O 



), and consequently we may average the exponen- 



tial of Eq. (|3j) to get 

P(t + T) 
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=-M(T)p(t) = S(T,T)p(t) 



(4) 



The equation 

M(r)p = (1 - 
n > 1 



for statistical equilibrium thus becomes 
I- jTL)p, which has the solution [pH] for 



Pn = PO 
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The overall constant p is determined by J2n°=o Pn = 1- 

Let the state of an atom emerging from the cavity be char- 
acterized by a binary "spin" variable s = ±1 where +1 de- 
notes the excited level. In statistical equilibrium the prob- 
ability 'P(s) of finding the atom in the state s after the in- 
teraction is given by T'(s) = u T M s p where u is a vec- 
tor with all entries equal to one, representing the trace over 
the cavity. If the detection efficiency is not 100% these for- 
mulas are modified accordingly jl5[]. The average spin is 
\i = (s) = V(+)—V{—) (and the atomic inversion — /i/2). Its 
variance is a 2 = ((s — /i) 2 ) = 1 — /i 2 . From the spin average 
we may in statistical equilibrium determine the average occu- 
pation number of the cavity mode (n) = rib + (1 — m) l^ r f^- 
The joint probability Vk(s%, S2) for observing the states of 
two atoms, with k unobserved atoms between them, is given 
by 



V k ( Sl ,s 2 )=u T M s *S k M 
1 + 7TL 



(6) 



where the matrix S = S(t,T) is defined in Eq. (Q). 
From this we derive the fc-step correlation function 
{ss)k = £ SllS2 sis 2 V k ( Sl , s 2 ) = V k (+, +) + V k (- -) - 
V k (+, -) - V k (-, +)■ It is easy to show that V k (+, -) = 
V k {— , +) using the stationarity condition for p. The normal- 
ized atomic correlation is 7^ = ((ss) k — (s) 2 )/a 2 . 

The state of the cavity can be characterized by the aver- 
age (n) of the occupation number operator h and higher- 
order observables such as the Mandel quality factor Qf = 
((n — (n)) 2 ) I (n) — 1, which has the property of being neg- 
ative for sub-Poissonian (and non-classical) states of the field 
tflfil]. We may further characterize the state of the radiation 
field by correlations between cavity observables at different 
times. If for example the occupation number h is measured 
twice, with k unobserved atoms passing in between, we con- 
sider the correlation function 



{nn) k — u T nS k hp 



(7) 



As before we define the normalized correlation to be j[ — 
((nn) k -(n) 2 )/((n-(n)) 2 ). 

If k is sufficiently large we define the asymptotic correla- 
tion lengths £a and £f by 
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which isolates the fc-dependence. The i?-factor secures that 
the correlation lengths are measured in units of physical time. 
That £ = £a = £,f follows from the fact that the time evolu- 
tion is governed by the same matrix in both Eqs. (^|) and (Q). 
The correlation length is therefore a convenient probe of the 
dynamics of the micromaser system. 

The generation of a sequence of outgoing atoms is 
a Markov process defined by the matrix S = (1 + 
7TL)~ 1 A/ (r), and correlations over long times are governed 
by its eigenvalues A„. They can be shown to be real, non- 
degenerate and bounded by 1 = Ao > Ai > ... >0 Jl7|]. 
Denoting the corresponding right (left) eigenvectors by 
we can decompose S as 
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The contribution to 7^ or 7^ from the A = 1 component of 
S k cancels out and it is the next-to-leading eigenvalue Ai that 
determines the correlation length through i?£ = — 1/ In Ai. 

In order to determine £ we have used two methods, a Monte 
Carlo simulation of the dynamical process and a direct nu- 
merical calculation of A„. We consider the maser transition 
63^3/2 «-> 61e?5/2 of 85 Rb with the single-photon Rabi fre- 
quency fl = 44 kHz. We choose the rate R to be 50 atoms/s 
in order to be close to the experimental situation of M\. In the 
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MC simulation we used a sample of 10 6 atoms for 100 differ- 
ent values of r (the corresponding real experiment would take 
a little less than a month to perform non-stop with 100% mea- 
suring efficiency!). We have compared the MC-data and the 
theoretical average of the atomic inversion as well as the Man- 
del quality factor. The agreement between theory and the nu- 
merical experiment is excellent and constitutes a consistency 
check of our computational method. The correlation length is 
measured by a least-squares fit to the data for 1 < k < K, 
where K is determined by the onset of noise which almost 
universally happens when the correlation has fallen to about 
0.1%. Typical values of K range from 10 to 250, depending 
on whether the correlation quickly or slowly drops into the 
noise. This leads to quite different reliability levels for the 
extraction of the values of the correlation length. Our results 
are presented in Fig. |l|. The error bars on the MC data rep- 
resent the precision of the fit to the exponential form Eq. (|8|), 
but do not include the systematic errors due to the difficulty in 
reaching the asymptotic region. 
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FIG. 1. Comparison of theory (solid curve) and MC data (dots) 
for the correlation length R£. The dotted and dashed curves corre- 
spond to subleading eigenvalues (A2.3) of the matrix S. The param- 
eters are those of the experiment in 

The other numerical method we have used is to find the 
roots of the equation det[M — A(l + 7TX)] = 0, by truncating 
the matrix to finite dimension. These roots are the eigenval- 
ues of S = (1 + r yTL)~ 1 Af. Writing the above equation as 
an eigenvalue problem for the matrix M — aL with a = A7T, 
we recognize that it is a Jacobi matrix (with non-vanishing el- 
ements only on the main diagonal and its two neighbouring 
subdiagonals). The off-diagonal elements are strictly positive 
for non-negative real a. It is known JT^ ] that such a matrix 
(for fixed a) has only real and non-degenerate eigenvalues. 
We have shown that the eigenvalues of S also have these prop- 
erties, and that A„ < 1, which is required for the convergence 
of the stochastic process. Numerical calculations of A„ con- 
verge rapidly for matrices of dimension above 50. They also 
agree well with the MC simulation, as can be seen in Fig. |IJ 



For r ~ 20 — 30 fis the MC result predicts shorter correlation 
lengths than the theoretical calculation, but we can see that 
they agree with a subleading eigenvalue. We surmise that the 
reason for this is that contributions of the subleading compo- 
nents of S are large for those values of k that can be reached 
in the MC calculation. In order to extract the actual long-time 
correlation length in this parameter region, it would be neces- 
sary to go much larger k and to increase the statistics consid- 
erably. Similar difficulties are expected to be encountered in 
a real experiment. 

The example of the correlation lengths shown in Fig. [ljhave 
typical parameters of the damping and the flux corresponding 
to a real experiment The amount of structure found in this 
figure is at first glance rather surprising. The mean photon 
number shows comparatively little structure in this region. We 
shall now explain the origin of structure in this figure, and 
refer to Fig. g, when necessary, which depicts the logarithm 
of the correlation length for larger flux values. The natural 
scaling variable in the large flux limit is 6 — gr^/N which is 
the one used in Fig. 0. Details of the theoretical derivations 
will be presented in [O]. 
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FIG. 2. The logarithm of the correlation length (in units of 
I/7 = N/R) for rib = 0.15 as a function of 9 = gr\f~N for various 
values of N — 10, 20, . . . , 100. The vertical lines indicate where 
new local maxima in the photon distribution function appear. Notice 
that for 9 > 8\ ~ 4.494, where the photon distribution acquires two 
maxima, the logarithm of the correlation length grows linearly with 
N for large N. 

In the first part of Fig. [l], for < r < 20 (as, the correlation 
length shows a peak around r = 7.2 fis which corresponds 
to the maser transition at 9 = 1. For smaller values of r the 
field of the cavity has thermal statistics with average number 
of photons independent of the flux N(= R/j) in the large 
flux limit. Above the peak, the average photon number grows 
proportionally with the flux. This peak, which is also clearly 
present in /1 and Qf, is well described by a semiclassical ap- 
proximation [p|Jl~^]. In Fig. |2| we see that in this region the 
correlation length remains constant in the large flux limit ex- 
cept at the maser transition point where it actually grows like 
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VN. 

At r ~ 20 fis, there seems to be crossings of eigenvalues 
in Fig. |l[ Looking closer at these points, we find that there is 
no actual crossing, consistent with the result mentioned above 
that there are no degenerate eigenvalues. On the other hand, 
the corresponding eigenvectors, in this case and p^ 2 \ do 
cross in the sense that the eigenvector corresponding to A2 
after the crossing is very close to p^ before the crossing. This 
has no analogue in fi or Qf. From the MC simulation we also 
see that the system may remain dominated by a subleading 
eigenstate as r passes through the transition point. In Fig. ^ 
we see how this crossing point as the flux increases creeps 
closer to 9\ ~ 4.494, where the photon distribution function 
gets a second maximum [ |r7[ ] . 

Above t — 32 fis the correlation length increases dramati- 
cally. The reason for that can be traced back to the appearance 
of two local maxima in the photon distribution function, and 
the tunneling time between them, as indicated in [[J and cal- 
culated explicitly in jl7|]. We have verified that, as the flux R 
of atoms increases, the maximal correlation length grows ex- 
ponentially in this region. For a value of R — 200/s (which is 
merely 4 times the flux in the experiment of [^]), the correla- 
tion length at the leading peak extends to 300 times the decay 
time of the cavity {i.e. a whole minute!). Such a violent be- 
haviour is not reflected in fi or Qf [PJt9|]. For these quantities 
the most dramatic behaviour is a finite jump in fi and a nar- 
row peak of finite height in Qf at the second maser transition, 
which coincides with the exponential peak in £. 

The second maser transition occurs when there is a jump 
in the position of the global maximum of the photon distribu- 
tion function. For TV = 10 this happens at r ~ 51.7^s. The 
three high peaks at r around 40, 50 and 70 fis have, however, 
a different explanation. They are related to the well-known 
trapping states J20| ] which occur when q n (r) vanishes. This 
happens for r^/n = kit with integer k. For vanishing n/,, this 
effect truncates the equilibrium photon distribution Eq. (||) at 
a definite photon number. Finite n\, smoothes out this effect, 
but it is still amply visible in Fig. |jj The dominant trapping 
peaks are numerically given by r = 41.2 fis, 50.5/zs, 71.4/xs 
which agrees very well with Fig. |l[ In Fig. ||we see how these 
peaks are rapidly suppressed relative to the leading and expo- 
nentially growing tunneling peak for larger flux. The trapping 
peaks occur at constant positions in r but move away as y/N 
in 9, whereas the tunneling peak remains at fixed 9. 

In the semiclassical regime there is a clear connection be- 
tween the long correlation at the maser transition at 9 = 1 
and the large value of Q f, i.e. super-Poissonian photon statis- 
tics. For 6 > 9\ ~ 4.494, and for large flux, Q f likewise 
has peaks at the same positions as £, namely when there are 
jumps in the position of the global maximum of the photon 
distribution. In this sense there is a clear relation between £ 
and photon statistics in the cavity. On the other hand, at 8\ for 
large flux the Mandel factor remains negative while £ starts 
growing exponentially, and this growth reveals the appearance 
of the second local maximum in the photon distribution which 
is not reflected in the single-time observables. 
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